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It was reported that Bose-Einstein condensation induces a spontaneous magnetization in spinor 
bosons. This phenomenon is called Bose-Einstein ferromagnetism (BEF). We propose a phenomeno- 
logical model consistent with the prediction of BEF and show that BEF might be attributed to the 
intrinsic magnetic moment of particles. Taking BEF into account, the phase diagram of an optically 
trapped spinor condensate is reexamined. 
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Alkali-metal atoms, such as 87 Rb, 23 Na, 7 Li, in which 
Bose-Einstein condensation (BEC) has been experimen- 
tally realized^], have an internal degree of freedom at- 
tributed to the hyperfine spin F. Accordingly, the or- 
der parameter of Bose condensed alkali atoms could have 
2F + 1 components. Thus, although the superfluid prop- 
erties of this system are considered to be essentially the 
same as those of superfluid 4 He, one can expect that the 
hyperfine spin degrees of freedom could bring about re- 
markable differences between the two systems. 

Earlier experiments leading to BEC in alkali atoms 
were performed in magnetic traps. A relatively strong 
external magnetic field H was applied to confine the BEC 
system. Because the atomic spin direction adiabatically 
followed H, the spin degree of freedom was frozen. As a 
result, the alkali atoms behaved like scalar particles al- 
though they carried spinQ- In 1998, Stamper-Kurn et 
ai. succeeded in cooling 23 Na in a purely optical trap and 
achieved BEC|3j. With H sufficiently small, the spin de- 
gree of freedom can become active and the spinor nature 
of the condensate can be manifest ed^jj]. 

The optical trap provides great opportunities to inves- 
tigate spinor Bose gases experimentally. Furthermore, it 
stimulates enormous theoretical interest in studyingvar- 
ious spin-related properties, e.g. the magnetism^, u| Q, 
of the gases. For an F = 1 Bose gas, the effective inter- 
actions V(r) between the atoms were derived asf" 



V(r) = -IffifylM* + y^F^ • F ctV ^ 7 , (1) 

where ip<r{ r ) is the field annihilation operator for an atom 
in hyperfine state \F, a) (a = 1, 0, —I) at position r. Re- 
peated sub-indices represent summation taken over all 
the hyperfine states. Besides a spin- independent scat- 
tering, a Heisenberg-like exchange interaction between 
hyperfine spins appears in the Hamiltonian, where the 
components of P = (F x , F y ,F z ) are 3x3 spin matrices. 
According to the pioneering work of Ho[j| and Ohmi and 
Machida 6] , this term determines the phase diagram of 
the spinor condensate: the ground state can be either fer- 
romagnetic or "polar" , depending on whether C2 is neg- 
ative or positive, respectively. However, recently Eiscn- 
berg and Lieb (referred to as EL hereafter) claimed 



that this scenario ignored an intrinsic ferromagnetism of 
spinor bosons. 

The intrinsic ferromagnetism in spin-1 ideal bosons was 
predicted by Yamada 9] and recently studied in the con- 
text of atomic Bose-Einstein condensation by Simkin and 
Cohen[l3. They calculated the magnetization M(H) 
and found that once BEC takes place, M remains fi- 
nite even if H = 0, suggesting the system is magnetized 
spontaneously. This phenomenon is called Bose-Einstein 
ferromagnetism (BEF). Caramico DAuria et al. pointed 
out that the coexistence of BEC with BEF also appears 
in interacting bosons [ll|. Moreover, rigorous proofs pre- 
sented by SiitoO and EL@ show that a polarized state 
is among the degenerate ground state of spinor bosons. 
EL studied a Bose model with a spin-independent, totally 
symmetric interaction v between particles, 
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where denotes the spatial coordinate of particle i. The 
specific model of Caramico D'Auria et al. has the same 
symmetry as the above Hamiltonian. 

In this communication, we propose a phenomenologi- 
cal model for the spinor Bose condensate. First of all, 
this model shows in a simple way that BEF can really 
take place at the same critical temperature as BEC. The 
results obtained are consistent with available theories. 
Then the phase diagram of optically trapped alkali atoms 
is reexamined, taking BEF into account. Furthermore, 
we try to understand the nature of BEF by analyzing 
the phenomenological model. 

To begin with, we derive an appropriate Ginzburg- 
Landau free energy density functional for second-order 
phase transitions that describes the coexistence of BEC 
and BEF. Generally, such a free energy density consists 
of three different parts: ft = fb + f m + /c> corresponding 
to the Bose condensed phase, the ferromagnetic phase 
and the coupling between the two phases (we omit the 
free energy of the normal Bose gas). Following Ginzburg 
and Pitaevskii|l3|. the free energy density of an isotropic 
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spin-F Bose-Einstein condensate is given as 

/ 6 = ^V*t.v*- a |*| 2 + ^|*| 4 
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where ^ 



(* J )* = is the com- 

plex order parameter and I*! 2 = St't'S'. Here the or- 
der parameter is defined as ty a = (ipa), but we note 
that other forms of the order parameter are possible 
(see LeggettQ)- The free energy density should have 
the same symmetry as the underlying Hamiltonian. The 
first three terms describe a system without any spin- 
dependent interactions, as given by Eq. (2). The fourth 
term has SU(2) symmetry as does the c-i term in Eq. (1), 
and describes the spin-exchange interaction. 

The free energy density of the ferromagnetic phase can 
be expanded in powers of M 2 : f m = aM 2 /2 + bM 4 /4, 
where M is the order parameter of the ferromagnetic 
state and is proportional to the magnetization density. 
In the usual theory of phase transitions, a changes sign at 
the ferromagnetic transition temperature and b > 00. 
Here, we assume that M is induced by the BEC itself, 
rather than by any other additional mechanisms, so f m 
should have a minimum at M = 0. That is, the ferromag- 
netic state would disappear if it were not coupled to the 
Bose condensed state. Therefore, f m (M 2 ) > / m (0) = 
for all values of M, so both a and b are positive. Near 
the critical temperature only the leading term need be 
considered, and f m takes the form, 



fm = \M 2 



(4) 



Physically, f m should contain the energy of the magnetic 
field inside a magnetized medium We suppose that 
the condensate only couples linearly to the ferromagnetic 
order, 
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where g is the coupling constant. This is the simplest 
coupling term that satisfies the physical situation, as in- 
dicated later. Our discussions are confined to the ther- 
modynamic limit and to one domain. Equations (4) and 
(5) comprise the basic ansatz of this work. We assume a 
and g are independent of the temperature. 

The total free energy for the Bose condensate can be 
written as Ft = J ftdV. Taking the variation with re- 
spect to ^ a and M, we derive 
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Eq. (6a) consists of a set of equations with respect to the 
spin index a. Eq. (6b) ensures that the magnetization is 
induced by the condensate itself. 

Possible phase transitions of the spinor Bose gas can 
be investigated by solving Eq. (6). For a spatially uni- 
form system, the first term in Eq. (6a) can be dropped. 
We take the bosons to have spin F = 1. It has been 
suggested that the condensate should be ferromagnetic 
even without the spin-exchange interaction\H. Ifll llOL lllj. 
To treat this case, we first set (3 S = 0. Two nontrivial 
solutions deserve further discussion |l6j. 

(I) Normal Bose-Einstein condensates. This solution 
corresponds to BEC without spontaneous magnetization. 
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with the free energy density 
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In this case, we have |*i| 2 = l^-il 2 but cannot deter- 
mine the proportion of l^il 2 to l^ol 2 - Thus, the system 
might be in "polar", "planar" or "equal spin" states[r?j|. 

(II) Bose-Einstein ferromagnets. This solution de- 
scribes BEC with a spontaneous magnetization: 



l*il -JZZ ' l*o| =l*-i| =o, 



(9) 



In this case the Bose condensate is fully polarized and 
the free energy density is 
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Comparing Eqs. (10) and (8), it can be seen that 
the BEF solution is indeed energetically favored in spite 
of a positive magnetization energy f m . The magneti- 
zation results in a net energy decrease by // — f\ l ~ 
(jr 2 Q! 2 /(2a / 9 2 )[l8|, suggesting that M ^ occurs spon- 
taneously. This magnetization depends on the onset of 
BEC, with its magnitude proportional to the condensate 
density, as pre dicted in previous theories Hoi Hl| . 

As usual|l3|. a can be expanded in powers of T c — 
T near the Bose condensation temperature T c : a(T) = 
a'(T c — T); (3 and a 1 are positive constants. According 
to Eq. (9), 
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The results are consistent with those obtained from a 
microscopic interacting boson model0|, suggesting that 
the phcnomcnological model describes the behavior of the 
spinor condensates quite well. 
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The realization of BEC in optically trapped alkali- 
metal atoms makes the study on BEF not only theoreti- 
cally interesting, but also experimentally attainable. We 
now turn our attention to such a specific BEC system. 
In this case, the (3 S term should be turned on. Using the 
notation (¥i,¥ ,*-i) = (4>ie ie \ fo, 0-ie 19 - 1 ), the p. 
term becomes / ?■,</>£ [<£? + 4>-i + 20i0_i cos(#i + + 
Ps{4>i -4> 2 -i) 2 /^ .19]. Without considering BEF, only this 
term determines the relative phases of the three compo- 
nent condensate wave functions: 9\ + 6-\ = tt for (3 S > 
and 6\ + 9-\ = for f3 s < 0. It therefore determines 
the ground state of the condensate: it is ferromagnetic 
for f3 s < and "polar" for (3 S > 00. However, this 
phase diagram is modified by BEF. After some algebra, 
the free energy density reduces to ff — —a 2 / (2(3) for 
"polar" states and to ff = -a 2 /[2(j3 + (3 S -g 2 /a)] for 
ferromagnetic states. Consequently, the phase diagram 
illustrated in Fig. 1 results, with a phase boundary at a 
finite value of f3 s , given by 



a 



(12) 



Equation (12) implies that there is a critical (3 S value 
leading to the "polar" state. 

Strictly speaking, the Landau expansion holds mathe- 
matically only near T c . In order to get the ground state 
phase diagram, we need decrease the temperature to zero, 
when the free energy density can be directly derived from 
the Hamiltonian via the Gross-Pitaevskii approximation. 
Because of the existence of BEF, we argue that the en- 
ergy due to the magnetization, as expressed in Eqs. (4) 
and (5), should be included in the T = free energy den- 
sity functional. Then the Gross-Pitaevskii equation for 
the F = 1 bosons described by Eq. (1) reads: 
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Here fi is the chemical potential. Equations (13) have 
similar forms as the phenomenological equations (6) , but 
their parameters have different meanings_20j, so here a 
and g are marked with a star. From Eq. (13), a phase 
diagram similar to Fig. 1 can be produced and the 
phase boundary becomes (ci) c = (g*) 2 /a* . The Gross- 
Pitaevskii approximation is quite accurate in describing 
the ground state of Bose gases. 

We now discuss the origin of BEF. The possibility 
that BEF arises from an explicit spin-exchange inter- 
action, like the C2 term in Eq. (1), has been ruled 
out@, 0, 0|. Our phenomenological model shows 
clues that BEF may be attributed to intrinsic magnetic 
moments of particles. Since spinor bosons carry mag- 
netic moments m, an effective magnetic field B = /xqM 
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FIG. 1: Phase diagram of spin-1 Bose condensates. The 
schematic /3 s -dependence of the two order parameters, "J/] 2 
and M, is outlined. Both decrease with increasing (5 a in the 
ferromagnetic state. At the phase boundary, |\&| 2 is contin- 
uous but M drops suddenly to zero, signaling a first-order 
phase transition. The intermediate region denotes the ferro- 
magnetic state attributed purely to BEF. 



is created inside a polarized body, owing to the superpo- 
sition of the intrinsic field of all the magnetic moments. 
fio is the permeability of vacuum and M = m| , I'i| 2 is the 
magnetization. Naturally, f m = £> 2 /(2/io) = HoM 2 /2 
can be regarded as the energy of the effective field. On 
the other hand, the magnetic moment of the condensed 
bosons does respond to the magnetic field regardless 
if the field is external or internal, so the spin direc- 
tion should follow the internal field B. Eq. (5) just 
is the general form of the coupling between the con- 
densate density and the internal magnetic field: f c = 
-^ "iM|^i| 2 = -M-B = -/i A/ 2 . Within this inter- 
pretation, f m + f c = — HqM 2 /2, which is just the self- 
energy of a magnetized body with magnetization M|l5|. 
Contrarily, in the "polar" or "equal spin" states the mag- 
netic moments of the particles are fully compensated: 
M = 0. Then at the mean-field level, an internal refer- 
ence particle cannot sense the moments of other particles, 
and the magnetic self-energy is zero. Obviously the fer- 
romagnetic state is energetically favorable at T = 0. At 
finite temperature, the ferromagnetic state can be de- 
stroyed due to an entropy increase. However, once Bose- 
Einstein condensation takes place, the temperature inside 
a condensate is absolutely zero and the ferromagnetic 
state remains stable. This accounts for the dependence 
of BEF on the onset of BEC. We note that such a inter- 
pretation coincides with the prediction of Yamada and 
Caramico D'Auria et al, because in their calculations of 
magnetization at a given field, the magnetic moment of 
spinor bosons was included. 

Simkin and Cohen attributed BEF to the Bose- 
Einstein statistics ^(J- However, we notice that although 
the Bose-Einstein statistics does not exclude a ferromag- 
netic state in spinor bosons, it cannot lead to a spon- 
taneous polarization. According to EL, the "fully" po- 
larized state is among the degenerate ground states but 
the degeneracy cannot be lifted spontaneously if we do 
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not consider the magnetic moment of particles. Rojo[2l| 
studied a two-component bosonic system which can be 
mapped into a "pseudo-spinor" Bose gas with F = 1/2 
(without magnetic moments) and showed that the ground 
state is not polarized when the interaction is component- 
independent. 

Within the interpretation that BEF is due to the mag- 
netic moment, the magnetic dipole interaction between 
particles is implicitly involved. The dipole interaction 
is caused by the intrinsic field of magnetic moments 
which forms the internal magnetic field when magnetic 
moments are arranged parallel. The influence of mag- 
netic dipole interactions on the properties of BEC has at- 
tracted considerable interest recently [l^- Although very 
weak, they can modify the ground state and the collective 
excitations of trapped condensates significantly 0, H^- 
Our approach gives a new perspective for studying the 
effect of the dipole interactions on the boson gas. It is 
worth noting that dipole interactions can hardly result 
in a ferromagnetic state in a real Fermi gas, because the 
interaction is too weak to approach the Stoner threshold. 

The intermediate area in Fig. 1 denotes the "pure" 
Bose-Einstein ferromagnetic state. This provides a direct 
way of detecting BEF experimentally. Once we have one 
experimental example pertaining to this regime, i.e., C2 > 
and the ground state is still ferromagnetic, the existence 
of BEF could be confirmed. The parameter ci cx (0,2 — 
ao)/3 can be estimated from experiments 5] , where a-i 
and ao are s-wave scattering lengths corresponding to the 
total spin two and zero channel, respectively. According 
to the above arguments, 

(cb) c = Mom 2 , (14) 

It is a rather small value. This makes the direct exami- 
nation of BEF a difficult task. 

Up to now, all of our discussions are devoted to charge- 
less bosons. One can expect that charged bosons exhibit 
more complex properties. For example, there would be a 
competition between the Bose-Einstein ferromagnetism 
and the Meissner effect, even though there is no exter- 
nal magnetic field applied. The latter tends to reduce 
or even to kill the magnetization. So, the nature of the 
ground state for charged spinor bosons is in our view still 
an open question. Relevant topics have been studied in 
triplet superconductors [23 ■ 

In conclusion, we have proposed a phenomenologi- 
cal model consistent with the prediction that magnetic 
bosons undergo a spontaneous magnetization together 
with Bose-Einstein condensation. We remark that this 
phenomenon deserves special attention due to two rea- 
sons: (i) It is an intrinsic property of spinor bosons and 
may be attributed to the magnetic moment of particles, 
(ii) The spontaneous magnetization is parasitic to the 
Bose-Einstein condensation, reflecting that the entropy 
of the Bose condensate is zero. 
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